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A measurement of branching fractions of the D° and D + mesons into the 4> meson is reported. 
The result is based on a data sample of 22.3 pb _1 collected at the CM energy of 4.03 GeV with the 
BES detector operated at the BEPC e + e~ storage ring. From tagged DD pair events the average 
branching fraction for a mixture of D° and D + is determined to be _B(D — + <f>X) = (1.29 ± 0.51 ± 
0.12)%. Upper limits at 90% confidence level are set to be B(D° -> </>X) < 2.5%, S(D+ -> <j>X) < 
5.0%, and B(D+ -» <j>e + v) < 1.6%. 



I. INTRODUCTION 



In an era of high precision experiments such as the B factories and the LHC, accurate measurements of b-flavored 
particles can benefit from a better knowledge of charm decays and their branching fractions. The inclusive decay 
D — > has not been measured Q. This branching fraction can serve as an independent check of the existence of 
additional exclusive decays of D mesons that contain a <fi meson [[[] , and for B° physics studies that use the (j>l pair to 
tag the Bg meson In addition, this branching fraction would be helpful in understanding the charm meson decay 
mechanisms. 



In this paper, we report a first measurement of the inclusive <f> decay branching fractions of charged and neutral D 
mesons and a new search for the exclusive semileptonic decay D + — > <f>e + v. 



> 

in 
o 

Cn ', II. DATA SAMPLE AND ANALYSIS METHODS 

o 

■ This measurement is based on 22.3 pb _1 of data collected in e + e~ annihilations at \fs = 4.03 GeV at the BEPC 
. during the 1992-1994. The BES detector has been described in detail elsewhere ||. 

D \ At y/s—4.03 GeV charm mesons D° and D + are produced via 

e+e" — > D+D~, D°E>o, 
CD. D+D*-,D*+D-,D°D*o 
T ^~i \ D* + D* _ D*°D*° 

> ■ 

followed by cascade decays of the D* mesons. However, the D* _ can decay either to 7r~D° or 7r°(7)D~, so that 

■ reconstructing a D meson does not necessarily determine whether the recoiling D meson is charged or neutral. In 
order to measure specifically B(D° — > (f>X) and B(D + — > <^X), the numbers of neutral and charged D mesons recoiling 
against a reconstructed D meson, and the type of the D meson from which the <p mesons come, must be determined. 
To this end two methods have been developed and are used to measure the inclusive branching fractions of the D 
mesons. 

A. The D° and D + combinative double tag method (CDTM) 

To measure inclusive <f> branching fractions of the D° and D + mesons, the is searched in the recoil side against a 
fully reconstructed D meson, and the numbers of events against the D° and D + decays, N* , , are determined, 

tag tag 

which can be related via 



throughout this paper, charge conjugation is implied. 
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N* =eN^o B{D~ -» 0X) + eNgo B(D°->0X), (1) 

tag ta S ta S 

N^ + = eNp+ B(D~ — + 0X) + eNj^t B(D°— >0X), (2) 

tag ta S tag 

to the branching fractions of their decays, B(D~ — ► 0X) and B(D° — > </>X), where N° , Nj^o , NE* + , and NEJ+ are 

-'-'tag U tag D ta g D tag 

respectively the numbers of D~ and D° decays on the recoil against D + and D° tags, and e is the detection efficiency 
of the </>. The values of N°„ , N°„ , N°+ , and N° + are determined from a measurement of the total production 

tag tag tag tag 

cross-sections of ractions e+e~ -> D*D*, D*D at 4.03 GeV by BES §. 

B. The recoil charge method 

At a/s =4.03 GeV, D*D* and D*D pairs are produced with no additional charged tracks. Charged pions arising 
from direct D* decays are very slow, and are mostly undetected in the BES detector. As a result, only decay products 
of the D + and D° are visible for most events. Let Qd be the charm flavor of the reconstructed D meson, and Q re c 
be the total charge of tracks recoiling against this D meson. The Q> ec distribution for D° (D + ) centers at (1), and 
has a spread of ±1. The recoil charge method selects neutral and charged D mesons according to 

Q rec = 0, or Q rec = Qu = -1 for D° tags (3) 

and 

Qrec ■ Qd < for D + tags (4) 

For inclusive D decays, the efficiency and the misidentification rate are 0.74±0.02 and 0.25±0.02, respectively, as 
obtained from Monte Carlo simulations, and are approximately the same for both charged and neutral D mesons. 
These numbers are confirmed using kinematically selected data events e + c~ — * D + D~ and e + e~ — * D°D°. For events 
in which a D tag and a recoil <fi has been fully reconstructed, the efficiency of the recoil charge method is improved 
over that of the inclusive D events. A Monte Carlo study of various D decay modes into final states containing a 
<f> has been performed, and the variations among their efficiencies are included in the systematic errors. For these 
events, the recoil charge method selects D meson type correctly 0.91±0.01±0.02 of the time, and misidentifies a D for 
0.09±0.01±0.02 of the events, where the first error is due to Monte Carlo statistics, and the second is systematic. 

III. DATA ANALYSIS 
A. Reconstruction of D, <f> Mesons 

Charged tracks are required to have good helix fits which have a normalized chi-square of less than 9 per degree 
of freedom. These tracks must satisfy | cos6*| < 0.8, where 9 is the polar angle, and be consistent with coming from 
the primary event vertex. For charged particles, a particle identification procedure is applied. A combined particle 
confidence level calculated using the dE/dx and TOF measurements is required to be greater than 1% for the 7r 
hypothesis. For the kaon hypothesis, > L n , where L is the likelihood for a particle type, is required. 

Charged and neutral D mesons are reconstructed via decays D° — > K~7r + , K~ir~ir + ir + and D + — > K~7r + 7r + . To 
reduce combinatorial backgrounds, only D mesons from e + c~ — > DD*,D*D* reactions are selected with cuts on the 
momenta of Kn7r combinations. Figures 1(a), 1(b) and 1(c) show the invariant mass distributions for events that 
pass the selections. The signals are fitted, and after having accounted for double counting, the number of D events is 
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determined to be 9054 ± 309 ± 416, where the first error is statistical and the second systematic. These D events are 
used as tagged c+c~ — > DD*, D*D* events in which the recoil side contain an unbiased D decay. 




1.7 1.8 1.9 2 2.1 

FIG. 1. Invariant mass distributions for K~tt + (top), K~7r + 7r + (middle), and K~ % + tt + ir~ (bottom). 

Table 1 summarizes the numbers of neutral and charged D mesons in the recoil against the reconstructed D tags. 
The averages from the CDTM method and the recoil charge method, calculated assuming a full correlation between 
their statistical errors, are 6803±303±322 and 2251±77±112 for D° and D+, respectively. 

The <p meson is reconstructed through its decay to K + K~. Figure 2 shows the invariant mass distribution of K + K~ 
pairs selected. Using convoluted Breit-Wigner and Gaussian functions plus a third order polynomial background to 
fit the mass spectrum, a mass of 1.0194 ± 0.0002 GcV/c 2 and a total of 1108 ± 70 4> events are obtained. In this 
measurement, a (f> signal window is defined as the region from 1.00 to 1.04 GeV/c 2 , as indicated by the arrows in 
Figure 2. 



TABLE I. Numbers of neutral and charged D mesons on the recoil 



method 


number of 


number of 




D events 


D + events 


CDTM 


6839±308 


2215±70 


recoil charge 


6767±297 


2287±83 


Average 


6803±303 


2251±77 
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FIG. 2. Invariant mass distribution of K + K pairs. 
B. Inclusive D — ► 

Figures 3(a) and 3(b) show the invariant mass distributions of K + K~ pairs from D+ and D°, respectively, as 
identified by the recoil charge selection criteria. The Kn7r invariant masses for the single tag are within ±2.5ctm d 
of the D masses. In this measurement, K + K~ pairs with masses in the ranges 0.98 - 1.00 GeV/c 2 and 1.04 - 1.15 
GeV/c 2 are taken as background for the <j>. The Knn mass regions from 1.7 to 2.1 GeV/c 2 , excluding regions within 
±3<tm d 01 the fit D masses, are defined as background control regions for the D mesons. As shown in Figures 3(a) 
and 3(b), 15 events are found as D0 candidates, and 14 events are selected as background outside the </> mass region. 
Using the D sideband events, a total of 0.5±0.5 background events has been estimated as the background among the 
D candidates. Subtracting the background contributions to both the D and the <f), we obtain an excess of 10.2±4.0 
events in the <j> signal region. 

The two D type identification methods, CDTM and the recoil charge method, are applied to these events to extract 
the numbers of <j> from specific D° and D + decays. Subtracting backgrounds estimated using the <f> and D side bands, 
the two methods determine 3.7±4.7 (CDTM) and 9.7±4.2 (recoil charge) D° -» </X events, and 6.5±5.5 (CDTM), 
and 0.5±1.7 (recoil charge) D + — > <pX events, respectively. Averaging over the two methods and assuming a complete 
correlation in their statistical errors, the number of D° — ► <pX and D + — ► <j>X events are set to be 6.7±4.5, and 3.5±3.6, 
respectively, and are used to determine their branching fractions. 
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FIG. 3. K + K invariant mass distributions for D° — > (fiX (top), and D + — > (fiX (bottom). 



C. Search for the decay D — > (fie v 



Among the 15 <fi candidates observed in the recoil side of the events, 4 are accompanied by at least one charged 
track which are within \cos9\ < 0.85. Each of these tracks is checked for consistency with being an electron using the 
dE/dx information. This electron identification requires that electron confidence level to be greater than 1%, and 
L c > L^. None of the accompanying tracks is identified as an electron. 



IV. RESULTS 

Assuming 10.2±4.0 signal D — > (fiX. events, and correcting for cf> meson detection efficiency of 0.084±0.006 obtained 
from a Monte Carlo simulation, the average branching fraction for the BES mixture of D° and D + mesons is measured 
to be 



B(D -> 4>X) = (1.29 ± 0.51 ± 0.12)%, 

where the first errors are statistical and second systematic. 

Based on 6.7±4.5 D° — > <fiX and 3.5±3.6 D + — > events, as determined in the previous section, 90% C. L. upper 
limits are set on specific D°, D + decays to be 
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B(D° -» (f>X) < 2.5%, 



B(D+ -» 0X) < 5.0% 

The results include systematic errors arising from uncertainties (±0.05%, ±0.06% and ±0.04%) in the numbers of 
singly tagged D mesons due to the choice of a background function and fit interval for the single tag samples and 
uncertainties (±0.08%, ±0.13% and ±0.09%) in the inclusive <j> efficiency. The combined effect of these sources is 
obtained by adding the uncertainties in quadrature, which yields total systematic errors of ±0.10%, ±0.14% and 
±0.10% for the D°, D+, and their sum, respectively. 

Based on zero candidate D + — > <pe + v events, and a detection efficiency of 0.0652, a 90% C. L. limit is set for the 
decays at 



B(D+ -» (pe+X) < 1.6%. 
V. CONCLUSION 

In summary, the inclusive branching fractions of the charged and neutral D mesons into a <f> have been directly 
measured. Comparing with the sums of the existing measurements on the exclusive D° and D+ decays containing a 
4> in the final states, these BES branching fraction values indicate little room for additional 4> decay modes of D° and 
D + mesons. 
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